Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





XMM- Newton light curves of the low-mass X-ray binary 
EXO 0748-676: dips, eclipses, and bursts 

Jeroen Homan 1 , Rudy Wijnands 2 , and Maureen van den Berg 1 

1 IN AF/Osservatorio Astronomico di Brera, Via Bianchi 46, 23807 Merate (LC), Italy; homan@merate.mi.astro.it, 
vdberg@merate.mi.astro.it 

2 School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews Fife, KY16 9SS, 
Scotland, UK; radw@st-andrews.ac.uk 

Received / Accepted 

Abstract. We present an analysis of XMM-Newton light curves of the dipping, bursting, and eclipsing low-mass X- 
ray binary EXO 0748-676, focusing on the variability on time scales of seconds to hours. The observed variability 
can be roughly divided in three types: dips, eclipses and bursts. We find that the appearance of the latter two, 
depends strongly on the strength of the first. We show that the absorption dips change from spectrally hard to 
spectrally soft as they become deeper, which supports suggestions that the source is composed of a spectrally hard 
compact source and a spectrally soft extended source. The fast variability in the soft light curve indicates that 
the large structures causing the dips are made up of smaller absorption cores. We present the first clear detection 
with XMM-Newton of eclipses below 2 keV in this source, and show that dipping activity is apparently unrelated 
to the source luminosity. We also test several proposed models for the spectral evolution during dips and confirm 
the presence of a scatter/reflection component in the eclipse spectrum. 
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1. Introduction 

Among the ~150 known low-mass X-ray binaries (Liu 
et al. 2001), there are 11 whose light curves show irregular 
and sometimes periodic intensity dips. These dips usually 
are ascribed to obscuration of the central X-ray source 
by thickened regions in the outer disc, that result from 
the impact of the accretion flow from the companion star 
(White & Swank 1982; Walter et al. 1982). Four of the 
dip sources also show (partial) eclipses, which are due to 
obscuration of the central X-ray source and/or the accre- 
tion disc corona by the companion star. The observed dips 
and eclipses are related to the inclination at which those 
sources are seen, which is probably between 75° and 90° 
(Frank et al. 1987), depending on the exact properties of 
the two stars and the orbital parameters. Because of their 
inclination and the regular obscuration of the bright inner 
parts of the system, these sources provide a good oppor- 
tunity to study the vertical and large scale structure of 
accretion flows in low-mass X-ray binaries. 

The transient source EXO 0748-676 is currently the 
only active low-mass X-ray binary that shows both dips 
and eclipses. Although the source was discovered in 1985 
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during an EXOSAT slew (Parmar et al. 1985), it was de- 
tected in Einstein archival observations from May 1980 
(Parmar et al. 1986); the source has probably remained ac- 
tive ever since, with 1-20 keV fluxes between ~ 10~ 12 erg 
cm -2 s _1 (with Einstein; Parmar et al. 1986) and 1.5 10 -9 
erg cm~ 2 s~ x (peak of the 1985 outburst). Observations of 
type I X-ray bursts (Gottwald et al. 1986, 1987) demon- 
strate that the compact object is a neutron star. The 
eclipses and dips revealed an orbital period of 3.82 hr 
(Parmar et al. 1986). Depending on the stellar parameters, 
the eclipse duration (8.3 minutes) and the orbital period 
give a system inclination of 75-82°. With EXOSAT, dip- 
ping activity (up to 80% decrease in 2-10 keV intensity) 
was observed at orbital phases ^0.8-0.2 and ~0.65 (phase 
zero being the center of the eclipse), with shallower dips 
at other phases. ASCA observations revealed that at low 
energies (1-3 keV) dipping (up to 100%) is present around 
most of the orbital cycle (Church et al. 1998). 

Observations with the Rossi X-ray Timing Explorer 
(RXTE) revealed the presence of quasi-periodic oscilla- 
tions (QPOs) with frequencies of ~0. 4-3.0 Hz and 695 Hz 
(Homan et al. 1999; Homan & van der Klis 2000). The ~1 
Hz QPO was only observed at low luminosities; a similar 
type of QPO was also seen in two other high-inclination 
systems (Jonker et al. 1999, 2000). This QPO is probably 
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caused by an opaque structure on top of the accretion disc 
that orbits the central source at a distance of ^1000 km, 
or by variations in the inner coronal flow as hydrodynam- 
ical simulations of the accretion disc corona show (Melia 
et al. 1991). The 695 Hz QPO, observed during a short 
period of increased luminosity early 1996, had properties 
similar to the kHz QPOs observed in many other neutron 
star sources (van der Klis 2000). 

Using the Reflection Grating Spectrometer (RGS) on 
board XMM-Newton, Cottam et al. (2001) recently found 
emission and absorption features whose properties also 
suggest the presence of a thickened accretion disc with 
plasma orbiting high above the binary plane. In their 
0.3-2.5 keV light curves no eclipses were seen at the 
expected times, which was confirmed by the apparent 
lack of eclipses in the 0.5-2 keV band of the XMM- 
Newton EPIC-MOS and EPIC-pn observations analyzed 
by Bonnct-Bidaud et al. (2001). Based on their results, 
the latter authors proposed a model in which the source 
is a superposition of a compact (~ 2 10 8 cm) spectrally 
hard Comptonizing corona and a more extended (~ 3 10 10 
cm) spectrally soft thermal halo. Spectral variations are 
then caused by variations in the absorption that affects 
the hard spectral component more than the extended soft 
component. This is in contrast with previous models of 
the source consisting of a compact black body and ex- 
tended Comptonized component (Church et al. 1998) or 
an absorbed plus un-absorbed Comptonized component 
(Parmar et al. 1986). 

In this paper we present an analysis of XMM-Newton 
light curves, focusing on variations on time scales of sec- 
onds to hours, both at high (>2 keV) and low (< 2 keV) 
energies. 

2. Observations and analysis 

For our analysis we used all XMM-Newton observations 
of EXO 0748-676 that were publicly available at the time 
of writing. Although XMM-Newton (Jansen et al. 2001) 
can observe simultaneously with different instruments, in 
many of our observations (all performed during the cal- 
ibration and performance verification phase) data from 
one or more instruments are absent or incomplete. We de- 
cided to use either EPIC-MOS1 (Turner et al. 2001) or 
EPIC-pn (Struder et al. 2001) data, depending on which 
of these cameras had the longest exposure. A log of the ob- 
servations, including the instruments from which the data 
were used, is given in Table 1. Note that the instruments 
were used in a variety of modes and with different filters. 
An analysis of observations 1 and 2 has previously been 
published by Bonnet-Bidaud et al. (2001). 

The MOS and PN data were reprocessed following the 
procedures described in the XMM-Newton data-analysis 
guides 1 and using the latest calibration data. Data were 
extracted from circles centered around the source, with 



Table 1. Log of XMM-Newton observations and data of EXO 
0748-676 used in this paper. From left to right: observation 
number for this paper, observation date with start and end 
time, the camera (EPIC-MOS1 or EPIC-pn) with mode and 
filter, and exposure time. 



Obs. 


Date & Start-End Time 


Camera 


Exp. 




(TT) 




[Mode,Filter] a 


(ks) 


1 


2000-04-04 17:19- 


-22:53 


MOS1 [FF,T] 


20.0 


2 


2000-04-21 03:59- 


-09:04 


PN [FF,M] 


18.1 




2000-04-21 09:50- 


-10:44 


PN [LW,M] 


3.2 




2000-04-21 15:22- 


-20:24 


PN [SW,M] 


18.1 


3 


2001-02-03 17:32- 


-19:32 


MOS1 [FF,M] 


7.2 




2001-02-04 01:23- 


-11:54 


MOS1 [FF,M] 


21.6 


4 


2001-10-13 00:19- 


-01:01" 


MOS1 [FF,M] 


32.9 



a Mode: FF = Full Frame, LW = Large Window, SW = Small 
Window; Filter: T = Thin, M = Medium 
b Observation ended on 2001-10-14 



radii of 90" (MOS1) or 40" (pn, in order to avoid the 
edges of the CCDs). Light curves were produced with a 
time resolution of 1 s and 10 s in the 0.3-10, 0.3-2 and 
2-10 keV energy bands. The 10 s light curves are shown in 
Fig. 1. A colour was defined as the ratio of the count rates 
in the 2 10 keV and 0.3 10 keV bands, and colour curves 
were created by taking the running average of 30 s inter- 
vals. Background light curves were extracted in a similar 
fashion (same radius, same CCD), with the center of the 
extraction circle at a location were no source contribution 
is expected - background extraction from an annulus cen- 
tered on the source, as is common practice, could not be 
performed, since the source was often located too close to 
the edge of the CCD. Variability in the background was 
never strong enough to warrant the exclusion of part of 
our data (see Fig. 1, gray lines). Note that the count rates 
shown in Figs. 1-4 are not background subtracted. 

Due to the brightness of the source, significant pile-up 
occurred in our observations, especially in the last obser- 
vation, where a clear hole could be seen in the center of the 
point spread function. To investigate the effects of the pile- 
up on the shape of light and colour curves, we made light 
curves by extracting data only from an annulus around 
the source position which excluded the piled-up core of 
the point spread function. We found that after scaling to 
the original count rates, within the statistics, those light 
and colour curves were indistinguishable from those made 
using the uncorrected data. In view of our statistics, we 
decided to use the initial light curves and not correct for 
pile-up. Because of the pile-up effects, which have a more 
serious effect on spectra than on light curves, we decided 
not to perform a detailed spectral analysis, except for the 
last the exposure of obs. 2. Details of this analysis, which 
was performed with the aim of testing the spectral models 
of Church et al. (1998) and Bonnct-Bidaud et al. (2001), 
can be found in the appendix. 



1 The XMM-Newton SOC User's Guide of the XMM- 
Newton Science Analysts System and An Introduction to 



XMM-Newton Data Analysts by the NASA/GSFC XMM- 
Newton GOF 



Observation 1 - 2000 April 4 



Observation 2 - 2000 April 21 
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Fig. 1. XMM-Newton light and colour curves of EXO 0748-676 for the 2000 April 4 (a) and 2000 April 21 (b) observations. Each panel shows roughly one orbital phase, 
with, from top to bottom, the colour curve, the 0.3-2 keV light curve, and the 2-10 keV light curve. For reference we also plot the background light curves (gray lines). The 
light curves have a time resolution of 10 s, whereas the colour curves show a 30 s running average. For clarity, error bars on the curves were omitted; they are in the order w 
of V0.1 x count rate. Type I X-ray bursts are indicated by Roman numerals. 



Observation 3 - 2001 February 3/4 



Observation 4 - 2001 October 13/14 



^ 1 

_o 
o 

<-> 

m 10 

m TJ 5 




-i i i i i 


MT - 


^ 1 

o 

o 

o 






i i i i i i i i i 


Soft 
cts/s 








tail 




- i i— 


r 


i i i i i i i i 


10 

P\ 
o o 

I-S 5 







0.5 0.6 


0.7 0.8 0.9 


1.0 1.1 1.2 1.3 1.4 1.5 




0.5 0.6 0.7 0.8 


0.9 


1.0 




1.1 1.2 1.3 


1.4 1.5 



<-> 

,n 10 



=_ 
10 - 

l-g 5 - 



w 



H h 



H 1 1 h 



H h 



2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 



<3 o = i i i 



H h 



H 1 1 h 



"i-g 5 - 




H h 



H h 



10 - 
°l 5 




3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 



1 



o 

<~> = I I I 



H h 



H h 



^ » 10 

o 5 

10 

X71 5 



H \ 



o t== 




4.5 4.6 4.7 4.1 



4.9 5.0 5.1 5.2 5.3 5.4 5.5 

Phase 




10 
5 

10 

5 






1.5 1.6 1.7 



.9 2.0 2.1 2.2 2.3 2.4 2.5 




2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 







l» ; 




Ixn 

s 

— I I I I I - 



1.6 6.7 6.8 6.9 7.0 7.1 
Phase 



7.2 7.3 7.4 7.5 



Fig. 1. . . .continued. Light and colour curves for the 2001 Feb 3/4 (c) and 2001 October 13/14 (d) observations. The sharp dips in the first orbits of panels c (just before 
the eclipse ingress) and d (around phase 1.05), during which no counts were observed, are probably instrumental effects. 
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Fig. 2. Folded 0.3-2 keV light curves of observation 2 (a) and 
observation 3 (b). The data were rebinned to 1000 bins per 
orbital cycle and, unlike Fig. 1, the count rate is plotted loga- 
rithmically. Type I bursts are still included. Weak eclipses are 
visible around phase zero. 



3. Results 

In Figs, la-d we show the light curves of EXO 0748-676 
in the 0.3-2 keV (soft) and 2-10 keV (hard) energy bands, 
together with curves for the colour. Count rates are plot- 
ted versus orbital phase, with phase zero assigned to mid- 
eclipse. Note that, because of different instruments, filters, 
and CCDs, the count rates of the four observations should 
not be compared directly; the low count rates in observa- 
tion 1 are due to the use of the thick filter instead of the 
medium filter, whereas the high count rates of observa- 
tion 2 are due to the higher sensitivity of the EPIC-pn 
detectors. 

As was already noted by Bonnet-Bidaud et al. (2001), 
at first sight the soft and the hard light curves have little in 
common; the soft light curves show a lot of rapid activity, 
whereas the variations in the hard light curves, except 
for eclipses and type I X-ray bursts, are much slower and 
usually weaker. 

3.1. Eclipses 

While the RGS and EPIC data analyzed previously 
(Cottam et al. 2001; Bonnet-Bidaud et al. 2001) were char- 
acterized by the absence of eclipses at low (< 2 keV) en- 
ergies, our fourth observation (Fig. Id) is the first XMM- 
Newton light curve of EXO 0748-676 to show clear eclipses 
at low energies. In contrast to those in the soft band, 
the light curves in the hard band always show eclipses. 
It should be noted that in the folded light curves of obser- 
vations 2 and 3 (shown in Fig. 2) during the phases of the 





-0.02 -0.01 0.01 0.02 

Phase 

Fig. 3. An enlargement of part of the folded 0.3-10 keV light 
curve of observation 4. The time resolution is 2 s. The dashed 
line represent the average count rate from 200 data bins cen- 
tered around phase zero. The dotted lines represent the start 
of eclipse ingress and the end of eclipse egress. 



eclipse the count rates in the soft band are actually below 
the average count rates for these observations (but not as 
low as some of the dips). Although not as apparent as in 
observation 4, this probably indicates that in these obser- 
vations the effects of the eclipse are present. In both bands 
the eclipses are not total, with residual count rates of ~5- 
8% and ~4-5.5% in, respectively, the soft and hard band 
of the fourth observation. Typical ingress/egress times in 
observation 4 are ~8-10 s (see Fig. 3), both in the soft and 
hard band. The folded 0.3-10 keV light curve of observa- 
tion 4 (Fig. 3) also reveals the presence of structures at 
the end of the eclipse ingress and beginning if the eclipse 
egress, similar to the 'shoulders' reported by Parmar et al. 
(1985). 

3.2. Dips 

Like in the ASCA observations analyzed by Church et al. 
(1998), dipping extends around most of the orbital phase 
and is most clearly visible in the soft band - only obser- 
vation 4 shows several intervals in which the 0.3-2 keV 
light curve is dip-free. With dipping being present most 
of the time, the most eye-catching features in the soft 
band are periods during which the strength of the dip- 
ping is reduced. These decreases in the dipping are the 
soft flares reported by Cottam et al. (2001) and Bonnet- 
Bidaud et al. (2001). Periods of reduced dipping show 
two typical time scales; extended intervals, lasting typi- 
cally for several minutes to an hour (Fig. 1) and shorter 
changes during those intervals, with time scales of sec- 
onds (Fig. 4a). When present, the extended intervals of 
(reduced) dipping in the soft and hard bands are proba- 
bly related to the orbital phase; they appear at a relatively 
constant phase and are present for several consecutive or- 
bits. This is most clearly seen in observation 2. Note that 
the interval of reduced dipping around phase 0.25 in this 
observation seems to show a slow, possibly random, shift 
towards higher phases; taking the moment during ingress 
at which the 'saturation' level in the soft band is reached 
as a reference point, the ingress of the dipping interval at 
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50 100 150 200 



Time (s) 

Fig. 4. Light curves showing the fastest observed dipping 
events in the 0.3-2 keV (a) and 2-10 keV (b) bands. The 0.3-2 
keV data are taken from observation 2 around phase 2.2 and 
have a time resolution of 3 s; the 2-10 keV are taken from the 
same observation around phase 3.85 and have a time resolution 
of 5 s. 



phases near 0.3, seems to occur at a later (~3%) phase in 
subsequent orbits. 

The phase extent and strength of the (decreases in) 
dipping undergo changes between orbital cycles. Unlike 
the long intervals, smaller structures do not seem to recur 
in later orbits - at least not at a similar phase. The rapid 
variations in the count rate in the soft band can be intense, 
with the count rate changing by factors of more than 10 
on time scales of only a few seconds (see e.g. Fig. 4a) - 
they are most clearly observed in observations 2 and 4. 
Excluding type I X-ray bursts and eclipses, the dynam- 
ical range (defined as the maximum count rate divided 
by minimum count rate, within one observation) of the 
soft band is >110. Variations in the hard band, seem to 
be more modest; the dynamical range is only ^30, and 
strong/ violent short-time variations, like those seen in the 
soft band, are not observed. For comparison we show the 
fastest observed dips in the soft and hard bands in Fig. 4; 
with a dip ingress of more than 30 seconds the fastest dip 
in the hard band is considerably slower than those in the 
soft band, which can be as fast a few seconds. During the 
low count rate periods (including the eclipses) the source 
is detected well above the background in the 0.3-2 keV 
band, with minima (observations 3 and 4) of ~0. 1-0.3 cts 
s -1 , compared to background count rates of ~0. 01-0. 03 
cts s _1 . 

As expected, the hard band is less affected by the dip- 
ping than the soft band. Dips in the hard band are only 
(but not always) observed when the count rate in the soft 




0.5 1 1.5 2 2.5 0.2 0.4 0.6 0.8 




1 2 3 4 5 6 0.2 0.4 0.6 0.8 1 




0.3-2 keV Count rate (cts/s) Colour 



Fig. 5. Left panels: 2-10 keV count rate vs. 0.3-2 keV count 
rate, Right panels: 2-10 keV count rate vs. colour. From top 
to bottom: observation 1, 2, 3, and 4. Each point represents a 
50 s average. The two dashed lines in the left panels represent 
the average background count rates in the two energy bands. 
The eclipses can be seen as small detached blobs. Note that 
while the count rates in the soft band during the eclipse are 
similar to those during the dips, they are clearly lower for the 
hard band. 

band is at a minimum. They are mostly observed at the 
phases around 0.7 and 0.9, which are also the phases which 
show the least variability and the lowest count rates in the 
soft band. The relation between dipping in the soft and 
hard band is clarified when we plot the count rates in the 
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Hard/Soft Peak Ratio 

Fig. 6. Histogram showing the observed ratio of peak count 
rates of the hard and soft bands for the 12 observed type I 
X-ray bursts. See also Table 2 

Table 2. Peak count rates of the 12 observed type I X-ray 
bursts in the soft (0.3-2 keV) and hard (2-10 keV) X-ray 
bands, together with their hard/soft ratio. Indicated is whether 
the bursts occurred during intervals of reduced dipping or full 
dipping (corresponding to the lowest count rates in the soft 
band, and, as opposed to reduced dips, sometimes also visible 
as dips in the hard band. For bursts I, II, and VII this was not 
completely clear. See Figure 1 for the labeling (I-XII) of the 
bursts. 



Burst 


Pers." 


Soft peak 


Hard peak 


Ratio 


I 


R(?) 


9.8±1.4 


19.8±2.0 


2.0±0.4 


II 


R(?) 


8.8±1.3 


18.4±1.9 


2.1±0.4 


III 


F 


8.4±1.3 


42.0±2.9 


5.0±0.8 


IV 


R 


20.6±2.0 


43.0±2.9 


2.1±0.3 


V 


F 


5.0±1.0 


38.0±2.8 


7.6±1.6 


VI 


F 


8.6±1.3 


104±5 


12.1±1.9 


VII 


R(?) 


87.4±4.2 


162±6 


1.9±0.1 


VIII 


R 


56.0±3.3 


99.4±4.5 


1.8±0.1 


IX 


F 


8.0±1.3 


27.0±2.3 


3.4±0.6 


X 


F 


1.8±0.6 


12.4±1.6 


6.9±2.5 


XI 


R 


10.8±0.5 


24.4±2.2 


2.3±0.2 


XII 


R 


15.0±1.7 


30.4±2.5 


2.0±0.3 



a Type of persistent emission in the soft band during which 
the burst occurred; R=Reduced dipping, F=Full dipping 



hard and soft band against each other, as has been done 
in Fig. 5; high count rates in the soft band are only found 
when the count rate in the hard band is high (i.e. outside 
dips). Dips in the soft count rate are however not always 
visible as dips in the hard band. During dips in the hard 
band the count rate does not drop below the level found 
in eclipses, whereas in the soft band the source can occa- 
sionally be weaker than during eclipse (see also Fig. 2). In 
all four observations the spectrum initially hardens as the 
count rate decreases, until the hard band starts to be af- 
fected, after which the spectrum starts to soften (see also 
Fig. A. 2). It is also interesting to see that in observation 
4 the colour saturates at the highest count rates and that 
the spectrum during the eclipse is softer than this satu- 
ration level. The motion of the source along the tracks in 
Fig. 5 is smooth, in the sense that (except during eclipses) 
the source never jumps from one end of the track to the 
other. 



3.3. Type I X-ray bursts 

In total 12 type I X-ray bursts were observed, 8 of which 
occurred in pairs with a separation of ~10. 65-14. 65 min- 
utes. When comparing the X-ray bursts in the hard and 
soft band we find that some of them hardly show up in the 
soft band. To illustrate this, we give the ratios of the hard 
peak and soft peak count rates for all 12 bursts in Table 
2 and plot them in Figure 6. Count rates were obtained 
from 5 s running averages to reduce the scatter in the 1 
s light curves. As no correlation was found between total 
count rate and the ratio of the peak count rates in the 
two bands, we assume that the effects of pile-up on our 
results are negligible. The peaks in the hard and soft band 
did not occur simultaneously - the quoted count rates are 
those at the time of the hard peak. Taking the soft count 
rates from the peak in the soft band (which occurred be- 
tween 9 s before and 7 s after that in the hard band) results 
in a similar picture. A clear concentration (7 out of 12) 
of bursts is present around a hard/soft peak ratio of ~2, 
with the other bursts all having larger values. While due 
to the rapid nature of the dips, it is not possible to rule 
out that some of the burst occurred during short periods 
of reduced dipping (which have similar time scales as the 
bursts), it seems that the 5 bursts that have ratios larger 
than 3 all occur during dipping intervals. Note that the 
small peak around phase 1.9 in the hard band of obser- 
vation 1 was not included, since the nature of this peak 
was not clear. If this peak, which occurred during a dip in 
the hard band, also was a type I burst its hard/soft ratio 
would be even higher than that for burst VI. 

3.4. Fluxes 

In order to compare the fluxes of the four observations as 
well as possible, we fitted the spectra of the source dur- 
ing the parts of the observations with the highest count 
rates. The spectra were fitted in XSPEC with various mod- 
els. While most models provided good fits in the 5-10 keV 
range, we were not able to satisfactorily fit the low energy 
parts of the spectrum with a consistent model (see also 
Appendix). For all four observations the unabsorbed 5- 
10 keV flux was consistent with 9 10~ n ergcm~ 2 s _1 . Our 
fits suggest that differences at lower energies are likely due 
to the complex effects of absorption. 



4. Discussion 

While at first sight the count-rate variations at low and 
high energies in EXO 0748-676 seem to occur indepen- 
dently from each other (Fig. 1), they are in fact strongly 
related (Fig. 2) and likely due to the same mechanism. 
This is supported by the fact that the source moves 
smoothly along the track in Fig. 2, which suggests that 
most of the variations in the light curve are caused by 
changes in only one parameter. In this section we dis- 
cuss how changes in the absorption can account for the 
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observed variations as well as for the appearance of the 
eclipses and type I X-ray bursts. 

4.1. Dips, eclipses, and bursts 

Both Church et al. (1998) (ASCA observations) and 
Bonnet-Bidaud et al. (2001) (XMM-Newton observations) 
showed that the spectral evolution of EXO 0748-676 can 
be largely explained by the effects of variable absorp- 
tion progressively covering two (stable) emission compo- 
nents. In both models the extended component is spec- 
trally softer than the compact central source and as such 
both are consistent with the observed behavior. However, 
as we show in the Appendix, the real spectral composition 
is likely to be more complex than suggested by Church 
et al. (1998) or Bonnct-Bidaud et al. (2001). 

As absorption is intrinsically less effective at high en- 
ergies, dipping is most pronounced at low energies and 
therefore spectrally hard. This can clearly be seen in our 
observations 1-3 and in the ASCA light curves of Church 
et al. (1998), in which dipping in the low energy bands 
extends over a much larger phase than in the high energy 
bands; the intervals of rapid increases in the count rate 
in observations 1-3 and most of observation 4 represent 
periods in which the absorption is reduced by a signifi- 
cant fraction. One would expect the ingress of dips in the 
soft band to precede that in the hard band, an effect that 
is indeed seen in our light curves and even more clearly 
demonstrated in Fig. 5. We therefore prefer to refer to the 
intervals of increased count rate as periods of reduced dip- 
ping instead of flaring, since we think this gives a better 
description of what gives rise to the changes in the light 
curves. 

Bonnct-Bidaud et al. (2001) suggested that part of the 
variability in their 'flares' is due to intrinsic changes in the 
halo luminosity (possibly due to reprocessing of type I X- 
ray bursts). However, the fact that the count rates in the 
eclipse in observation 4 are similar to those in observation 
3, which had the same instrumental setup but was much 
less luminous outside the dips in the soft band, sugge sts 
(assuming that the same part of the halo was left uncov- 
ered during the eclipse) that the halo luminosity did in fact 
not vary significantly. Moreover, since the 'flares' occur at 
similar phases for several orbits, delayed reprocessing of 
type I X-ray bursts is basically ruled out as a trigger for 
halo-luminosity changes. 

The apparent lack of eclipses in the soft band of the 
first three observations can easily be explained; the addi- 
tional absorption column presented by the companion star 
has almost no influence below 2 keV, since most of the soft 
flux from the central source was already absorbed at the 
time of eclipse - the projected size of the halo also has to 
be significantly larger than that of the companion star to 
not be affected during an eclipse. The fact that the count 
rates outside of the eclipse are sometimes lower than those 
during the eclipse suggests that the absorbing structures, 
as seen from the central source, occasionally subtend a 



larger solid angle than the companion star. This indicates 
that in the vertical direction the absorbing material ex- 
tends well beyond the expected scale height of a thin disc. 

The behavior of the type I X-ray bursts can also be ex- 
plained in this framework; when they occur during intense 
dipping intervals, i.e. when absorption is high, the soft flux 
of the burst is more affected than the hard flux, resulting 
in ratios larger than 2. Those that occur during intervals 
of reduced dipping are affected less by absorption, result- 
ing in ratios of ~2, which is apparently the normal ratio 
for these bursts in this source. The spread in ratios larger 
than 2 probably reflects various stages of absorption. We 
note that many of the bursts analyzed by (Cottam et al. 
2002), which showed red-shifted absorption lines, occurred 
during periods of significant absorption below 2 keV. 

4.2. Location of the absorbing material 

Since we know that the absorbing material has a height 
of more than 8-15° above the orbital plane it has to be 
constantly refreshed in order to be visible at a more or 
less fixed position in the orbital frame, for at least several 
orbital cycles; if not, it would disappear on a time scale 
that is shorter than the orbital period, because the time 
scale for changes in the vertical scale height in the accre- 
tion disc is comparable to the dynamical time scale (Frank 
et al. 1992). Building on the work of Lubow & Shu (1976), 
Frank et al. (1987) proposed a model in which an accretion 
stream overshooting the point of impact with the disc set- 
tles in a ring-like structure at a radius of ~ 10 10 cm (for a 
1.4 M Q neutron star and P or b — 3.82 hr). As the result of 
an ionization instability, a two-phase medium is formed, 
which consists of a hot thin gas and cool clouds and is 
present between phases 0.3 (at which the material joins 
the disc) and 0.8. The presence of such a two-phase ab- 
sorber, consisting of both neutral and ionized gas, seems to 
be confirmed by recent Chandra observations of the source 
(Jimenez-Garate et al. 2003). If we apply the equations of 
Frank et al. (1987) to EXO 0748-676, we find a typical 
size for the cool clouds of 4 10 7 cm and a column density 
of 6 1 23 cm~ 2 . At a distance of 10 10 cm this should cause 
dips with a time scale of ~1 s, quite similar to the short- 
est observed time scales. While these clouds can explain 
the rapid variations in the soft band, they cannot (within 
this model) explain the occurrence of dipping at phases 
between 0.8 and 0.3. 

Another possible site of absorbing structures could be 
a thickened disc rim; Rctter et al. (2002) found that the 
dips in the high-inclination X-ray binary X1916-053 are 
responsible for the positive superhump period observed in 
the X-rays and optical. They suggest that a thickened disc 
rim can explain both the X-ray and optical superhump. 
Positive superhumps (quasi-periodicities a few percent dif- 
ferent from the orbital period) are observed in some cat- 
aclysmic variables (Warner 1995) and low-mass X-ray bi- 
naries (O'Donoghuc & Charles 1996) and are interpreted 
as a signature of an eccentric precessing accretion disc. 
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Evidence for such superhumps may be present in observa- 
tion 2 in the form of dipping intervals shifting in phase (see 
section 3.2). However, the phase shift of these dips might 
as well be a random effect; unfortunately, the data sets for 
the other observations either are too short or contain too 
many gaps to study this in more detail. 

Whatever the exact nature of the absorbing mate- 
rial is, the fourth observation shows that its structure 
and/or presence can change considerably on time scales 
shorter than eight months (which is the interval between 
observations 3 and 4). It is not clear what causes these 
changes. The 5-10 keV flux during the brightest intervals 
of each observation are consistent with each other, sug- 
gesting that, if proportional to the 5-10 keV flux, the mass 
accretion rate remained constant. A possible explanation 
could be the precession of a tilted accretion disc, evidence 
for which has been found in several other low-mass X-ray 
binaries (e.g. Larwood 1998, and references therein). A 
tilted accretion disk, and hence precession, may also be 
present in EXO 0748-676 as was suggested by Crampton 
et al. (1986) on the basis of optical observations. If true, 
the absorbing structures could remain unchanged, but de- 
pending on the phase of the precession more (or less) of 
the absorbing material moves through our line of sight. 

4.3. Constraints from time scales 

The ingress and egress time scales we observed during the 
eclipses are on the order of ~8-10 s, somewhat larger but 
consistent with values obtained with EXOSAT and RXTE 
(Parmar et al. 1986; Hertz et al. 1997). Depending on 
the density profile of the companion star's outer atmo- 
sphere (see a discussion on this topic in Parmar et al. 
1986), this means that 92-95% of the 0.3-10 keV flux and 
94.5-96% of the 2-10 keV flux originates in a region with 
a radius of 3.5 — 3.9 10 8 cm or smaller. The latter num- 
ber is obtained from the duration of the eclipse ingress 
and egress; depending on the properties of the secondary 
(see Parmar et al. 1986), binary parameters (such as sec- 
ondary's radius, binary separation and inclination) were 
derived, which, using the eclipse time scales, translated 
into upper limits on the size of the eclipsed central ob- 
ject. As was already found by Parmar et al. (1986), near 
the bottom level of the eclipse, the eclipse ingress and 
egress progress slower than at the higher count rate lev- 
els (Fig. 3). This cannot be due to the density profile of 
the secondary's atmosphere, and suggests a central source 
geometry that is at least partly extended. 

The shortest time scales in the soft band are in the 
order of a few seconds and those in the hard band are 
~30 s. These time scales are difficult to interpret, since 
the distance from the central source and velocity across 
our line of sight of the absorbing material are not well 
known. However, the time scales for the dips in the soft 
band, which are shorter than the eclipse ingress and egress 
times, suggest that the central region is smaller than the 
values given above. The shortest observed time scales for 



the two bands are quite different. This difference can be 
explained if we assume that the absorbing material con- 
sists of small clouds (patchy absorber) as in the model of 
Frank et al. (1987) discussed above, each with a density 
high enough to considerably affect the soft flux, but not 
the hard flux. As the absorbing material (containing the 
small clouds) starts moving into our line of sight, the soft 
band is immediately affected, but it takes more time to 
have the required number of clouds present in our line of 
sight to affect the hard flux. An additional effect of such 
a patchy absorber is that the variations in the soft band 
are much stronger than in the hard band - quite the op- 
posite of what is found for other types of variability in 
X-ray binaries, which usually get stronger towards higher 
energies. 

4.4. Residual emission 

From Figs. 1 and 5 it is clear that eclipses and dips in 
the soft and hard band both tend to saturate at a level 
that is significantly above the background. As discussed 
above, part of this residual flux can be contributed to 
(fluorescent) emission from the extended thermal halo 
(Cottam et al. 2001; Bonnet-Bidaud et al. 2001). However, 
as suggested by observations of type I X-ray bursts during 
eclipses in EXO 0748-676 (Parmar et al. 1985; Gottwald 
et al. 1987), -4% of the 2-6 keV flux from the central 
source is scattered into our line of sight during eclipse. 
This scattering occurs in the accretion disc corona, which 
might be the same component as the extended thermal 
halo, and accounts for most of the residual flux in the hard 
band; residual emission in the hard band during eclipses 
in observation 4 is also ~4%. The slightly higher fraction 
of residual flux in the soft band (compared to the hard 
band) can be contributed to additional emission from the 
extended thermal plasma present above the disc. 

5. Summary 

We have analyzed XMM-Newton light curves of the low- 
mass X-ray binary EXO0748-676. We find that the ap- 
pearance of the light curves in the soft and hard band can 
be well explained by progressive absorption of two spectral 
components. From the colour of the residual emission dur- 
ing eclipses, and from the colour evolution with increasing 
degree of absorption, we conclude that the source contains 
an extended source of X-ray emission that is spectrally 
softer than the compact source, as suggested earlier by 
Church et al. (1998) and Bonnet-Bidaud et al. (2001). 

In the soft band, the count rate during dipping inter- 
vals sometimes decreases below the level during eclipses, 
while the source is always detected well above the back- 
ground, both in the soft and hard band. This implies that 
(at certain phases) the absorbing material as seen from the 
central source, subtends a larger solid angle than the pro- 
jection of the companion star, but that it is smaller than 
the extended soft component. The recurrence of dipping 
structures in the light curves in consecutive orbits indi- 
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cates that the absorbing material can maintain a more or 
less fixed position in the corotating frame for at least sev- 
eral orbits. However, the drastic change of appearance in 
the light curve of October 2001 with respect to that of 
February 2001, indicates that the position and/or density 
of the absorbing material can vary within months. The 
difference in the shortest time scales for dipping activ- 
ity in the soft and hard band, suggests that the dipping 
is caused by a patchy absorber that affects soft energies 
more efficiently 

Future observations with XMM-Newton, in particu- 
lar for several continuous orbital periods, may shed more 
light on the possible presence of superhumps and on the 
evolution and structure of the absorbing material. 
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Appendix A: Spectral analysis 

Several attempts have been made in the past to model the 
spectral evolution of EXO 0748-676 as a function of in- 
creasing absorption. In this appendix we test the models 
by Church et al. (1998) and Bonnet-Bidaud et al. (2001) 
that were already discussed extensively in the main text. 
Although the soft-band light curves presented earlier in 
this paper and the recent spectral results by Jimenez- 
Garate et al. (2003) strongly suggests that the obscuring 
medium is far from being uniform, these models treat the 
absorber as having a smooth density profile. Moreover, the 
extended components that are part of the two models are 
assumed to have a uniform brightness. Notwithstanding 
these simplifications we feel that a vigorous test of the two 
models may yield important information on the structure 
of the central source and the absorber. 

A.l. Data reduction 

For our spectral analysis we used the small-window mode 
EPIC-pn data of obs. 2 (i.e. the third data segment), the 
only part of our data set that was not affected by pile-up. 
Source data were extracted from a circle with a 38" ra- 
dius, centered around the source; background data from a 
similar circle in the opposite corner of the window. Only 
events with pattern 0-4 were selected and data from bursts 
were discarded. Data were selected either on time or the 
0.3-10 kcV count rate; the times and count rates for the 
various selections can be found and seen in Table A.l and 
Fig. A.l. Note that our time selections include those made 
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Table A.l. Time and count rate selections used for our spec- 
tral analysis. These selections are also displayed in Fig. A.l. 
Note that the eclipse (E) is excluded from the LI count rate 
selection. 



Selection 


Time" 


0.3-10 keV count rate 




(TT) 


(cts/s) 


E 


16:58:40-17:05:20 


— 


Fl 


17:42:30-18:30:00 




Q 


18:30:00-18:48:40 




F2 


18:48:40-19:23:40 




D 


19:23:40-19:45:20 




LI 




1.0-6.0 


L2 




6.0-12.5 


L3 




12.5-17.5 


L4 




22.5-27.5 


L5 




25.0-45.0 



a Time on April 21 2000. 



by Bonnet-Bidaud et al. (2001) (Fl, Q, F2, D) with the 
addition of an eclipse spectrum (E). The count rate selec- 
tions, which do not cover the entire count rate range (see 
gaps in Fig. A.l) were converted to time intervals using 
a 10s 0.3-10 keV light curve. The resulting spectra were 
fitted with the Xspec package (Arnaud 1996, VI 1.2) with 
the models described below, using a canned response ma- 
trix. The ancillary response file was created using arf gen. 

Before we discuss our spectral fits, some comments 
should be made on the two selection methods. First, se- 
lecting spectra based on temporal behavior (e.g. dipping 
activity), as was done by Bonnet-Bidaud et al. (2001) (see 
selections Fl and F2), will result in mixing of spectra 
with rather different properties; as is clearly illustrated by 
Fig. 5, spectral proprieties can change considerably with 
count rate. Selecting (narrow) count rate intervals, as was 
done by Church et al. (1998), will reduce these problems 
a bit, but mixing still occurs in those intervals where the 
count rate does not correlate with colour. Second, from 
other low-mass X-ray binaries we know that spectral evo- 
lution of the central source can occur on time scales that 
are considerably shorter than the exposure we analyzed 
(i.e. a few hours). Spectra that were taken more than a 
few hours apart, and which therefore have likely differ- 
ent central source properties, could end up having similar 
count rates because of the varying absorption, resulting in 
the mixing of different spectral properties. Unfortunately, 
tracing changes in the central source spectrum would re- 
quire splitting our exposure in shorter pieces, resulting in 
a low statistical quality of the spectra. 

A. 2. Spectral fits 

The eclipse spectrum and the count rate selected spectra 
are shown in Fig. A. 2. For a display of the time selected 
spectra we refer to Fig. 4 of Bonnet-Bidaud et al. (2001). 
Fig. A. 2 clearly shows that, when going to lower count rate 
selections (from L5 to L3) the spectrum above ~3 keV re- 
mains largely unchanged, explaining the hardening of the 
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Fig. A.l. Light curve of the third data segment of obs. 2 (0.3- 
10 keV) showing the time and count rate selections applied for 
our spectral analysis. Type I X-ray bursts are not shown. The 
labels L1-L5 indicate the count rate level selections (see main 
text). The other labels indicate the time selections and stand 
for eclipse (E), flaring interval 1 (Fl), quiescence (Q), flaring 
interval 2 (F2), and dip (D). Note that the eclipse was not 
included in the LI count rate selection. 




0.2 0-5 12 5 10 
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Fig. A. 2. Energy spectra for the 5 count rate selections and 
the eclipse interval. From top to bottom: the spectra for the 
selections L5 to LI, and the eclipse spectrum (big filled square 
symbols). All spectra are rebinned for display purposes. 



spectrum (see Fig. 5). From L3 to LI that part of the spec- 
trum also starts to become affected, which, because of the 
saturation at low energies (compare L2 and LI spectra), 
leads to a softening of the spectrum. The eclipse spectrum 
nicely shows that the soft part of the spectrum is consis- 
tent with that of the deepest dips, whereas at high energies 
the source is almost completely obscured (compare eclipse 
and LI spectra). 

Spectral fits were made with the models used by 
Church et al. (1998) and Bonnet-Bidaud et al. (2001). The 
first model was only used in conjunction with the count- 
rate selected spectra (L1-L5), whereas the Bonnet-Bidaud 
et al. (2001) model was used for both types of selection 
(to study the effect of selection method). 
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Table A. 2. Spectral fit results. Models and energy ranges are given for each of the fits. Lower limits on ABB are 90% confidence. 
All other errors represent la confidence intervals. 



Model: Church et al. (1998) Model: Church ct al. (1998) + Gaussian 

Range: 0.8-10 keV xLj/d-O-f-=L01/1594 Range: 0.15-12 keV Xred/d-Q-f .=1.18/1916 





L5 


L3 


LI 


L5 


L3 


LI 


ABB ( a ) 


0.01(19) 


0.2(2) 


>32 


0.01(2) 


4.2(4) 


>91 


kT bb (keV) 


1.6(1) 






1.46(4) 








6.5(4)- 10 -5 






6.8(2)- 10 -5 






AGA ( a ) 


0.010(4) 






0.031(2) 






APL (°) 


0.4(4) 


3.2(1) 


14.0(4) 


0.12(1) 


2.43(5) 


13.7(2) 


f 


n nq+ ' 46 


0.90(2) 


0.955(2) 


0.4(3) 


0.853(1) 


0.963(1) 




1.33(1) 






1.31(1) 






N Pi n 


0.0181(4) 






0.0177(1) 






E H „e (keV) 








0.524(2) 






a Une (keV) 
N Kne ( d ) 








0.103(2) 












0.032(1) 







10 atoms cm' 
photons keV" 1 > 



-1 . b 



L39 / (Dio) , with L39 the luminosity (10 ergs s ) and Dio the source distance (10 kpc) 
at 1 keV ; d total photons cm -2 s _1 in the line 



A. 2.1. Church et al. (1998) model 

The Church et al. (1998) model consists of a black body 
component from the neutron star, plus Comptonized 
emission from an extended and progressively covered 
accretion disc corona (ADC): ABB*BB + APL*f*PL + 
AGA*(l-f )*PL, where BB and PL are the black body and 
power-law terms, AGA the Galactic absorption term, ABB 
and APL the variable black body and power-law absorp- 
tion terms, and f the progressive covering fraction of the 
ADC. Note that in this definition of the model ABB and 
APL include Galactic absorption and that the two PL com- 
ponents represent the covered and uncovered parts of the 
same spectral component. In all fits with this model the 
columns densities (ABB, AGA and APL) have a lower limit of 
10 20 atoms/cm 2 , which is a factor of ~10 lower than the 
average galactic column density in the direction of EXO 
0748-676 (Dickey & Lockman 1990). AGA and the param- 
eters for the black body and the power law were treated 
as single parameters for all spectra. 

We started by performing fits to the LI, L3 and L5 
spectra in the energy range used by Church et al. (1998) 
(0.8-10 keV). Fig. A. 3a shows the results of this fit; 
the best-fit parameters can be found in Table A. 2. The 
X 2 ed /d.o.f of 1.01/1594 is good, although ABB for spectrum 
L3 seems to be too low, given the covering fraction of 90% 
that we find. More worrisome are the unacceptably large 
residuals below 1 keV and also around 1.6 keV in the fit 
to spectrum LI. This part of the spectrum is not shown in 
Fig. 3 of Church et al. (1998), so we cannot check whether 
such residuals were also present in their fits to the ASCA 
data. Most likely, these residuals are due to the soft ex- 
cess found by Thomas ct al. (1997) and also discussed by 
Church et al. (1998). The presence of the this soft compo- 
nent is much clearer in the 0.15-12 keV spectra and a fit 
with the model in this energy range is clearly worse than in 
the 0.8-10 keV band (x 2 ed =1.76). Following Church et al. 



(1998) we add a Gaussian line to the model, which, like 
the power law component, is progressively covered. The 
best fit with this model is shown in Fig. A. 3b; the best-fit 
parameters can be found in Table A. 2. Again, even though 
the Xred = l-18 is acceptable, large (localized) residuals are 
present - this time in the fits to all three spectra. As a 
general remark we note, that since most data bins are in 
the rather featureless part of the spectrum above a few 
keV, which is fitted well by the models, the residuals be- 
low 2 keV (which is the part of the spectrum we are most 
interested in) have only a limited effect on the X 2 ed - 

When plotting the separate spectral components it be- 
comes clear that rather than being a narrow line, the 
Gaussian is very broad and fits a large part of the contin- 
uum below 1 keV. Moreover, the ADC component dom- 
inates the spectrum over the entire energy range in all 
three spectra. The covering fraction and the two separate 
absorption components (AGA and APL) act in such a way 
on the (very hard) ADC component that the uncovered 
power law (together with the Gaussian) fits the low en- 
ergy part of the spectrum, whereas the covered power law 
fits the high energy part. In the L5 spectrum the black 
body component only contributes 20% of the absorbed 
0.15-12 keV flux. 



A. 2. 2. Bonnet-Bidaud et al. (2001) model 

In the Bonnet-Bidaud et al. (2001) model the source con- 
sists of a spectrally hard compact central component and a 
spectrally soft extended component. These are modeled, 
respectively, by a power law (PL) a thermal Raymond- 
Smith model (RS), each affected by their own absorption 
component: APL*PL + ARS*RS, where APL and ARS both 
include the galactic absorption and are not allowed to be- 
come lower than 10 20 atoms/cm 2 . 
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Model: Church et al. (1998) Model: Church et al. (1998) + Line 

Spectra: L1, L3, L5 - Range: 0.8-10 keV Spectra: L1, L3, L5 - Range: 0.15-12 keV 




12 5 0.2 0.5 1 2 5 

Energy (keV) Energy (keV) 




Fig. A. 3. Model fits to various XMM-Newton spectra of EXO 0748-676. Spectra, models, and energy ranges are given above 
each panel. All spectra were rebinned for display purposes. In the panels a, b, and d spectra L5, L3 and LI are plotted from 
plot to bottom. In panel c the order of the spectra, from top to bottom, is F2, Fl, Q and D. See Table A. 3 for fit parameters. 



The first fits with the Bonnet-Bidaud et al. (2001) 
model were made to the time selected 0.15-12.0 keV spec- 
tra. The element abundances were fixed to the values 
found by Bonnet-Bidaud et al. (2001), while the other 
parameters were allowed to vary. The power law index, 
the power law normalization and the temperature of the 
thermal component were treated as single parameters in 
the joint fit to all four spectra. We obtain very simi- 
lar values as Bonnet-Bidaud et al. (2001), but with a 
slightly higher Xred °f 1-34 (d.o.f.=2904) compared to 
their Xr e( j/d-0.f.=1.27/2467. This difference can be partly 
attributed to a different grouping of the data (minimum of 
25 counts per bin instead of 20), higher maximum energy 
(12 keV instead of 10 keV), and also due to improved cal- 
ibration of the instruments. In any case, as is clear from 
Fig. A. 3a, large localized residuals and broad deviations 
are present; e.g. in the fit to the dip spectrum around 
2 keV, and in the fits to the flare and quiescent spectra 
above 3 keV. 

Bonnet-Bidaud et al. (2001) do not argue why they 
did not treat the normalization of the thermal component 



as a single parameter for all four spectra, as they did for 
the normalization of the power law component. It is not 
clear to us why the properties of an extended component 
would vary on a shorter time scale than those of a com- 
pact component. Their argument that the increases in the 
normalization of the thermal component are due to repro- 
cessing of energy from type I X-ray bursts is not supported 
by our observations of bursts and flaring intervals. The in- 
creases in the normalization occur with decreases in the 
absorption of the thermal component, suggesting that the 
variable normalization tries to compensate for the absence 
of a progressive covering factor in the model. Indeed, ty- 
ing the value of the normalization to a single parameter 
results in a xled of 1.92. 

Next we fitted the Bonnet-Bidaud et al. (2001) model 
(i.e. without tying the normalization of the thermal com- 
ponent, to stay consistent with the previous fit) to our 
intensity selected spectra LI, L3, and L5 with the abun- 
dances of N, O, Ne, and Mg left free to vary. Using the 
fit-results from the time selected spectra as a starting 
point we obtained a Xred ~l-6 with similar fit parame- 
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Table A. 3. Spectral fit results. Models and energy ranges are given for each of the fits. All errors represent la confidence 
intervals. 



Model: Bonnet-Bidaud et al. (2001) Model: Bonnet-Bidaud et al. (2001) 

Range: 0.15-12 keV Xred/d-Q-f -1.35/2904 Range: 0.15-12 keV xLj/d-Q-f-=l-54/1914 





Fl 


F2 


Q 


D 


L5 


L3 


LI 


apl n 


1.82(2) 


1.28(2) 


5.7(1) 


11.6(2) 


1.34(6) 


2.77(3) 


17.9(3) 




1.29(1) 








1.44(1) 








2.19(1)-10~ 2 








2.85(2)-10~ 2 






ARS ( a ) 


0.048(1) 


0.052(1) 


0.18(1) 


0.24(1) 


0.033(2) 


0.010(1) 


0.110(4) 


kTfls (keV) 


0.50(1) 








1.42(3) 








3.35(4)-10~ 2 


4.86(5)-10~ 2 


2.0(1)-10 -2 


1.98(6)-10~ 2 


3.79(1)-10 -2 


5.6(1)-10 -3 


2.77(5)-10~ 3 



a 10 22 atoms cm -1 ; b L39 / (D10) 2 , with L39 the luminosity (10 39 ergs s _1 ) and D10 the source distance (10 kpc) 

c photons keV" 1 cm~ 2 s _1 at 1 keV ; d W~ 14 / (4ivD'\) J n^nndV , with where Da is the angular size distance to the source 

(cm), rie is the electron density (cm -3 ), and uh is the hydrogen density (cm -3 ) 



ters. However, a lower x 2 re d °f 1-53 was obtained with a 
higher temperature and a much higher abundance of N 
(~80 times solar compared to ~10 times); the fit results 
can be found in Table A. 3 and are shown in Fig. A. 3d. 
Again large residuals are present around 2-3 keV in the 
LI spectrum. Unlike the fit to the time selected spectra, 
here the normalization and the absorption of the extended 
component do not correlate. 

A. 2. 3. Eclipse spectrum 

A powerful test of both spectral models can be made with 
the spectrum obtained during the eclipse. As the compact 
component is eclipsed, one only expects to see contribu- 
tion from the extended component. In the Church et al. 
(1998) model this is the (progressively covered) power law 
+ Gaussian combination and in the Bonnet-Bidaud et al. 
(2001) model this is the thermal halo. A fit with a power 
law + Gaussian at 0.52 keV gives xj! d /d.o.f.=1.5/18, 
whereas the fit with the Raymond-Smith model gives 
X 2 ed /d.o.f.=1.9/15. It should be noted that for the latter 
model we obtain both very high temperatures and abun- 
dances. Fixing these parameters to more reasonable val- 
ues (like those obtained from the fits to the non-eclipse 
spectra) gives a xfed values of 3.5-5, with large positive 
residuals above 2 keV. 

Given the strong indications that the scat- 
tered/reflected light from compact source/neutron 
star does contribute to the eclipse in this source - as 
suggested by observations of type I X-ray bursts during 
eclipse (Parmar et al. 1985; Gottwald et al. 1987) - it 
is not surprising that both models fail to fit the eclipse 
spectrum. We therefore include the compact components 
again in both models, keeping in mind that 1) their con- 
tribution should be around 5% of that measured outside 
eclipse (i.e. the percentage of scattered flux measured for 
the type I X-ray bursts in eclipse) and 2) that the Njj 
for this component should not be larger than that during 
dips, as the scattered light should be much less affected by 



the absorbing structures. Adding the black body (with kT 
fixed to 1.46) to the power law + Gaussian combination 
docs not lead to a significant improvement of the fit. 
However, adding a power law (with index fixed to 1.35 
and abundance all to 1) to the Raymond-Smith model 
leads to a much improved fit: x^ ed /d.o.f.=l. 12/19, with 
kT~0.4. The normalization of the power law component 
in this fit is ^2-3% of that in the fits to the non-eclipse 
data and the Nh is ~ 8 • 10 21 atoms cm -2 , in reasonable 
agreement with our expectations and also similar to the 
results obtained from eclipse spectra in MXB 1659-298 
(Sidoli et al. 2001). Assuming that a considerable part 
of the reflecting/scattering medium is obscured during 
the eclipse, it is quite likely that the contribution to 
the spectrum of this reflected/scattered component is 
considerably higher outside of eclipse, partly explaining 
the problems encountered when fitting those spectra. 

A. 3. Concluding remarks 

Our fits clearly show that current models are not able to fit 
the absorption affected spectra of EXO 0748-676 during 
dips. While the fits to the eclipse spectrum favor the ther- 
mal halo model of Bonnet-Bidaud et al. (2001) this model 
has severe problems with the deep dip spectra. Our eclipse 
spectrum provided strong evidence for the presence of an 
additional scatter/reflection component, which should be 
an integral part in future models. Such models should also 
include more advanced models for the absorbing material 
and also for the structure of the thermal halo. 



